Logarithmic diffusion and porous media equations: a unified description 
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In this work we present the logarithmic diffusion equation as a hmit case when the index that 
characterizes a nonhnear Fokker-Planck equation, in its diffusive term, goes to zero. A linear drift 
and a source term are considered in this equation. Its solution has a lorentzian form, consequently 
this equation characterizes a super diffusion like a Levy kind. In addition is obtained an equation that 
unifies the porous media and the logarithmic diffusion equations, including a generalized diffusion 
equation in fractal dimension. This unification is performed in the nonextensive thermostatistics 
context and increases the possibilities about the description of anomalous diffusive processes. 
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INTRODUCTION 



Diffusion is a very usual phenomena in the nature and 
it happens, in general, when the system goes to an equi- 
librium state . Therefore it is of fundamental relevancy 
in physics, chemical, and biological processes. The lin- 
ear dependence in the time growth of the mean square 
displacement {x'^{t)) oc t, or alternatively in the variance 
(when {x) 0), is the finger point of the Brownian move- 
ment and usual diffusion. It is a direct consequence of 
the central limit theorem and the Markovian nature of 
the underlying stochastic process In contrast, the 
anomalous diffusion is in general characterized by a non 
linear variance growth in the time, that is to say, the dif- 
fusion will be considered anomalous when the behavior 
goes out the former. This kind of diffusion has a fun- 
damental role in the analysis of a large class of systems 
such as plasma diffusion 0|, diffusion in turbulent fluids 
, fluids transportation in porous media |^ , chaotic 
dynamics 6] , non-gaussian behavior of the heartbeat , 
diffusion on fractals p, |^ 0, 0|, anomalous diffusion 
at liquid surfaces |l2l. in the study of vibrational en- 
ergy in proteins |l3l |. among other physical systems. In 
the anomalous behavior description the variance growth 
can be a power law kind ((Aa;)^) oc i'', or even present 
another pattern. In this classification, when 77 > 1, we 
have a superdiffusive process, ry < 1, a subdiffusive one 
and 77 = 1 describes usual diffusion. Furthermore, in an 
anomalous diffusive process, the variance can be not fi- 
nite, describing a Levy process, although it presents a 
well defined index that characterizes such process 
The description of this kind of processes is based in the 
validity of the generalized central limit theorem, termed 
Levy-Gnedenko, which states that, by iV-fold convolu- 
tion, a distribution with divergent lower moments tends 
to one of the Levy stable class . 

It is possible to simulate the anomalous behavior of the 
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diffusion applying generalizations on the ordinary diffu- 
sion equation. It can be performed by introducing an ap- 
Topriate time-dependence or spatial dependence 

l 0, in the equation's coefficients. Also we can 
apply fractional derivatives Jl, 22, 23] . However, the 
introduction of nonlinearities reveals a large set of pos- 
sibilities to describe anomalous diffusive processes. An 
interesting characteristic of the nonlinear Fokker-Planck 
equation is that its stationary solutions, and some par- 
ticular time-dependent solutions, are such that maximize 
the Tsallis entropy, a nonextensive entropic form pro- 
posed in the last years by Tsallis [23 . l25j| . This effort 
becomes necessary when the Boltzmann-Gibbs statistics 
fails, for instance in the presence of long range interac- 
tions or memory effects and fractal phase space structure 
(see Ref. |2^ for recent review). The nonextensive sta- 
tistical mechanics has shown a very fruitful scenario to 
study anomalous diffusion. In this way, several works 
dealing anomalous diffusion were developed in such con- 
text. The connection between this formalism and the 
nonlinear Fokker-Planck equation was first pointed by 
Plastino and Plastino 1271 . The results were enlarged by 
Tsallis and Bukman ^S\, including a linear drift term. 
A phenomenological microscopic dynamics of the nonlin- 
ear Fokker-Planch equation is presented in 29], as well a 
nonlinear Fokker-Planch equation with state-dependent 
diffusion^]. A Tsallis MaxEnt solutions of the nonlinear 
Fokker-Planch equation applied to the study of correlated 
anomalous diffusion |3lL l32 ] , nonlinear fractional deriva- 
tive Fokker-Planch like equation ^22, 2^ , aging in in non- 
linear diffusion (33| , anomalous diffusion with absorption 
[3^ , anomalous diffusion in a fractal dimension j 
are some examples among other references. In this direc- 
tion, the generalized thermostatistic, based on the nonex- 
tensive Tsallis entropy, becomes a natural scenario as the 
correlated anomalous diffusion as the anomalous diffu- 
sion like Levy (variance not finite) |35l Is^ . IstL l3^ . In this 
context, the index of nonlinearity joined to the nonlinear 
Fokker-Planck equation, in the correlated anomalous dif- 
fusion case, can be connected to the entropic index q of 
the Tsallis entropy and with the respective generalized 
distributions. This relation is just q — 2 — v, where u is 
the nonlinear index, and it is easy to verify that when 
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q = 2 the diffusive term becomes trivial. We will show 
that this difficulty can be overcame by a logarithmic dif- 
fusion equation, or, in another words, this equation rep- 
resents an alternative for the limit case where the expo- 
nent that characterizes the nonlincarity of the diffusive 
term goes to zero. In this way, the logarithmic diffusion 
equation is just inspired in nonlinear diffusion equations 
like Fokker-Planck 19, 27, 28] . We wiU show that the 
non-stationary solution for the logarithmic equation is a 
lorentzian, and this implies that the second moment is 
not finite. It indicates that this diffusion equation is re- 
lated to supcrdiffusive processes, of the Levy kind. On 
the other hand, in the nonextensive statistical mechan- 
ics scenario is possible to unify the correlated anomalous 
diffusion equation, or porous media equation, and the 
logarithmic diffusion equation. Thus we enlarge the de- 
scription of diffusive processes by a unified equation that 
congregates the correlated anomalous diffiision and a like 
Levy one. 

This work is structured as follow. In the section II 
we will present the logarithmic diffusion equation and its 
exact solution. In the section III this equation is solved 
including a linear drift term. In the following section 
it is solved when a absorption term is present. We will 
obtain in the section V the stationary solution, and in the 
section VI we will perform the unification of the porous 
media and logarithmic diffusion equations. It is included 
the connection with the diffusion on fractals. Finally, in 
the last section will be presented the conclusion and final 
remarks. 



II. LOGARITHMIC DIFFUSION EQUATION 

To motivate our discussion, we will first consider the 
porous media equation. 



dp 
dt 



It has been employed in the analysis of percolation of 
gases through porous media {v > 2) j^^i thin saturated 
regions in porous media {v = 2) |4?il |. a standard solid-on- 
solid model for surface growth {v = 3), thin liquid films 
spreading under gravity {v = 4) among other s [43| . 
A solution of Eq. ^ is the generalized g-gaussian |23| , 
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Z{t) 



[l-(l-g)/3(t)r 



21 — 



(2) 



By direct substitution of into ^ it is easy to verify 
that 



dt 
1 dZ 

'Z~dt 



2DI3vZ'i-'^ 



(3) 
(4) 



and we obtain (3{t) oc ^-2/(3-9) and Z{t) oc i^/^^-q) Qf 
course, when g — > 1 we recover the results for the usual 



diffusion. This solution is obtained taking account the 
relation q — 2 — u. And what happens when q = 21 It 
is obvious that in this case the right side of the equation 
vanishes all the time and ^ is not the solution now. 
However, Eq. |(5J, for q = 2, recovers a very important 
function: the lorentzian. A question arises: how would 
be the shape of the nonlinear diffusion equation whose 
solution was a lorentzian? 

It is known that In a; decrease more slowly than any 
power a;*", when r goes to zero positively, so we are influ- 
enced to substitute by In p in Eq. ^ when v ^ 0. In 
this perspective, we introduce the logarithmic diffusion 
equation 



dp ^ j^ d"^ Inp 
dt dx^ 



(5) 



This equation emerges in plasma physics |43| and, in par- 
ticular, has been predicted for cross-field convective dif- 
fusion of plasma including mirror effects ■43| . The same 
equation describes the expansion of a thermalized elec- 
tron cloud 45] and also arises in studies of the central 
limit approximation to Carleman's model of the Boltz- 
mann equation 46, 47]. 

It is verified directly that Eq. ((Sj presents the 
lorentzian solution 



p{x,t) = 



1 



1 



Z{t) [l + /3(i)a;2 



(6) 



The substitution of © in © shows that and Z{t) 
obey the equations 



1 dZ 

'z^'dt 



1 dZ 

'Z^'dt 

1 dl3 

]3Z~dt 



2DP 
-2DP. 



(7) 

(8) 



(1) They can be decoupled taking account the relation 



tions are 



1/2 



which is valid all the time. The solu- 



and 



l3it)=Ml + 2DPoZoty 



Z{t) = Zoil + 2Dl3oZot). 



(9) 



(10) 



It is interesting to point that for the lorentzian the vari- 
ance is not finite. This fact, characteristic of the Levy 
distributions, indicates that the logarithmic equation is 
associated to supcrdiffusive regimes. In fact, a dimen- 
sional analysis of the logarithmic diffusion equation (jsj 
and the Levy diffusion equation dp/dt = Dd^p/d\x\^, 
with p = 1, leads to the common ballistic behavior in 
the sense that x scales as t. By other side, Eq. |(SJ is a 
particular p = I case of 



dp d _^ dp 
dt dx ^ dx 



0<p<2. (11) 
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When the coefficient diffusion D{p) (or the thermal co- 
efficient of conductivity) can be approximated as p~^, 
its divergence for small p causes a much faster spread of 
mass (or heat) than in the linear case {p = 0), justify- 
in^ the terminology supcrfast diffusion to these processes 

For the A^-dimensional case, Eq. jSJ assumes the shape 



(12) 



with — X^i^i /dxj and whose lorentzian solution, 
p(|x|, t), presents 

m = /3o[l - 2{N - 2)Df3oZo i]^/^^-^) (13) 

and 

Zit) = Zo[l - 2{N - 2)Dl3oZo tf^^^-^\ (14) 

III. PRESENCE OF EXTERNAL FORCES 

Now we analyze Eq. © describing a process that in- 
clude a linear drift term (external force) F{x), 

For the case F{x) — ki — k2X, the solution of Eq. H15|l is 
the drifted lorentzian 

1 1 



p{x,t) = 



(16) 



^ ki - k2Xo 

= -4:D(3Z + 2k2 

^ 2D(3Z-k2. 



(17) 
(18) 



In fact, (3{t), Z{t), and xo(t) obey the equations: 

dxo 
~dt 

1 dZ 

^^-2DpZ-k2. (19) 

Equation H17|) does not depend of the index related to 
the non-linearity, therefore it arises not only in the usual 
diffusion equation but also in the porous media ones 7^ 
1). Thus, the solution of is 



Xo(t) 



xoiO) 



feat 



1 



-k2t 



(20) 



In this turn, the solutions of H18|l and H19|) have the form 



and 
where 

git)- 

with 5(0) = 



/3(t) = /3o [1 - 5W] 
Zit) = Zo[l- git)] , 



ik2 - 2Df^oZo) 



1) 



'k2t 



(21) 
(22) 

(23) 



IV. LOGARITHMIC EQUATION WITH 
SOURCE TERM 



We can consider the presence of a time-dependent 
source term. In this case, Eq. © is written as 



(24) 



The source term in this equation can be removed by an 
appropriate change in the solution 



pir,t) = e(- ''o p(r,i) 



(25) 



and we rewrite the time variable. Thus, the solution of 
(EH) is 



pix,t) 



ZiTit))[l + PiTit))x^ 



■ exp 



ait')dt' 



with /3(T(t)) and Zirit)) being of the form ® and ifTUIl . 
where D = 1 and t is replaced by rit) — Dit')dt' , 

with -D(t) = Dexp can to extend this 

solution for the iV-dimensional case by employing H13|l 
and ijnjl. 

V. STATIONARY CASE 

We define now the probability density current 

91np 



J = Fp-D- 



dx 



(26) 



(27) 



so that Eq. (|15|l represents a continuity equation: 

^ + ^ = 0. 
dt dx 

In the stationary case, dJ /dx = 0, and this implies 
J — cte — 0, since we are applying the condition that 
the current vanishes at infinity. This fact permits us to 
write, considering F = —dV/dx, 



dlnp 
dx 



1 dV 
13 dr' 



whose solution is 



l + f3V ' 



(28) 



(29) 



where it was assumed V^(0) = 0, /3 — po/ D and piO) — 
po = 1. Note that the structure H29|) is preserved in the 
A'^-dimensional case. 

The solution (|29|l is the q = 2 case of the generalized 
exponential 



p=[l-(l-g)/3y]i/(i-') 



(30) 



which, in this turn, recovers the Boltzmann distribution 
in the limit g — > 1. Furthermore, we remind that the 
Boltzmann distribution is the stationary solution of the 
usual diffusion equation, so H29|l is the analogous for the 
logarithmic equation. 
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VI. UNIFICATION OF POROUS MEDIA AND 
LOGARITHMIC EQUATIONS 

The unification of equations and H12(l will be accom- 
plished by using the generalized logarithmic function, the 
g-logarithm, defined as 



In, X 



ci-« - 1 



(31) 



In this way, the unified equation proposed in this work 
presents the following structure 



|=5V^ln 
at 



q-lP, 



(32) 



and it is clear that, when g — > 2, the function Ini recovers 
the logarithmic one. When q = 2 — u ^ 2, we have 
ln.q_i = {p'^ — l)/v and Eq. ^A2\ reobtains the porous 
media equation ^ with D = D/v. 

We can consider the non linear diffusion equation with 
radial symmetry, taking account the spatial-dependence 
in the diffusion coefficient, r~^, and a non integral di- 
mension d 18,] , 



with 



dr 



dr 



(33) 



(34) 



With the operator A written like this there are not re- 
strictions for the possible values for d. In this way, d can 
be interpreted as a fractal dimension in an embedding 
TV-dimensional space. Equation H33|) . in the v = 1 case, 
recovers the diffusion equation introduced in Refs. [IjQ- 
Analogously to (|32|l . we enlarge the applicability domain 
of Eq. 133|l and we have 



dp 
at 



DAln„ 



-1 



with A defined by 13411 . The ansatz 
1 



Z{t) 



\l-{l-q)mr 



A] — 



(35) 



(36) 



with A = 2-1-0, remains valid and substituted in Eq. H35|) 
conducts to the equations 



dZ{t) 

dt 

dp{t) 
dt 



DXd/3{t)Z'^{t) 
~DX^f3^{t)Z'^-\t). 



(37) 



Such equations are decoupled and solved and their solu- 
tions are 



(38) 



and 



with 



Z(i) = Zo[l + At]''/[-^+''(i-«)], 



D\[X + d{l-q)](3oZf 



' 



(39) 



(40) 



/3o — /3(0) and Zq = Z{0). Thus we obtain the stretched 
lorentzian if A < 2 and short if A > 2. In fact, for q = 2, 
the solution is 



1 



Zit) 



1 



1 -I- P{t)r^ 



(41) 



where (3(t) = (3o[l + At]-^^^^-"^^ and Z{t) = Zo[l + 
^^]d/[A-d\ Observe that when 6' = Oandrf = 7V = l 
we obtain the results l(2Jl for g 7^ 2, and Q for q — 2. 

We pointed that the factor (2 — q) implicated Eqs. ^ 
and (j^ became trivial at the value q — 2. With our 
strategy we eliminate this undesirable behavior. From 
the performed generalization in this section, we can con- 
jecture a nonlinear diffusion equation, that presents p'^, 
can be extended to = 0, substituting p'^ by lnq_ip, 
with q — 2 ~ ly. Thus, in a general way, we can unify all 
the equations presented in this work by the equation 



dp 
dt 



N 

E 



d 

dx. 



D,. 



dx 



■ ln„ 



-1 P 



N 



i=l 



(42) 



VII. CONCLUSION 

In this work we enlarged the application domain of 
the porous media equation for the v = case and we 
presented the logarithmic diffusion equation as an al- 
ternative procedure to this limit case. Its solution has 
a lorentzian form and it can characterize superdiffusive 
processes, of Levy kind. A unification of porous media 
and logarithmic diffusion equation is obtained, and in a 
more general form, with the fractal nonlinear diffusion 
equation. This unified description is performed in the 
nonextensive statistical mechanics scenario. This accom- 
plishment represents a progress in the formal description 
of diffusive processes and in its solutions as well. In this 
direction, the final proposed equation interpolates oth- 
ers ones with consecrated position in the literature. It 
is desirable that the equations presented in this work, 
or in particular special cases of them, show physics sit- 
uations in which there is competition among different 
mechanisms that generate anomalous diffusion. 
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